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Abstract 

Recent experimental studies strongly suggest the influence of glial purinergic 
transmission in the modulation of synaptic dynamics. By releasing adenosine 
triphosphate (ATP), which accumulates as adenosine, astrocytes tonically sup- 
pressed synaptic transmission. The delayed multi-step feedback of the glial 
adenosine with the neuron suggest the existence of hysteresis phenomena, which 
are investigated in the present study from the theoretical point of view. The 
model suggests that a memory operator, tripartite synaptic plasticity, governs 
the mysterious delayed feedback inhibition caused by the action of adenosine 
on neuronal Ai receptors and provides a powerful tool for further dynamical 
modeling tasks on tripartite synapses. 



1. Introduction 

The nervous system consists of two classes of cell, the neuron and glia. Ac- 
cording to the classical view of the nervous system, the numerically superior 
glial cells have inferior and passive roles in that they provide an ideal environ- 
ment for neuronal-cell function. Recent evidence has demonstrated, however, 
that there is a dynamic reciprocal feedback communication between glia and 
neurons at the synapses ([T], 0, [T^, [Mj)- Neurotransmitters released from 
presynaptic neurons evoke Ca^"*" concentration increases in perisynaptic glia. 
It is known that all types of glia (i.e. microglia, oligodendrocytes, astrocytes 
and Schwann cells) express a wide variety of neurotransmitter receptors f|18j). 
Activated glia, then release neuroactive substances, called gliotransmitters, in- 
cluding glutamate and ATP. These gliotransmitters feed back on the one hand 
onto the presynaptic terminal either to facilitate or to suppress further release 
of neurotransmitter, and on the other hand stimulate directly postsynaptic neu- 
rons by producing either excitatory or inhibitory responses ([1], [B], [7], [S], [10], 

ini, El). 

Released ATP from glia could excite neurons directly through activation of 
P2X receptors. Alternatively, released ATP, once converted to adenosine, could 
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inhibit neurons by activating Ai receptors, as it does in the retina. Adenosine 
could also act presynaptically via Ai and A2 receptors to either depress or 
potentiate synaptic transmission ([3], [5], [6], [9], [IT], [13], [14], [15], [16]). 

It has been observed that during hypoxia, astrocytes contribute to regulate 
the excitatory synaptic transmission through the release of adenosine, which 
by acting on Ai receptors reduces presynaptic neurotransmitter release f|13j). 
Because adenosine is implicated in the control of wake-to-sleep transitions as 
well as responses to hypoxia, the identification of the role of the glia in regulating 
this transmitter offers mechanistic insights into these processes ([2])- 

The kinetics of ATP hydrolysis and adenosine accumulation provides a frame- 
work with unique spatiotemporal conditions to control synaptic transmission. 
It takes « 200 ms before adenosine begins to accumulate. The delays (pre- 
accumulation delay, accumulation delay) and different states of the action of 
glial ATP on neurons, suggest further investigations to reveal the mysterious 
nature of the delayed multi-state process ([TT]). 

The aim of this article is the introduction of the concept of hysteresis in 
synaptic scale especially in the astrocytic ATP-related activities. This concept 
provide a simplified framework for the description of ATP kinetics that can 
be easily embedded into mathematical models dealing with neurochemistry of 
tripartite synapses. Thus, one can consider the present work as a tool for further 
investigations. 

Based on Preisach model for hysteresis, nonideal relays — Raj,bj are in- 
troduced. By /ij-weighted superposition of a finite bundle of relays the operator 
over the half plane P of admissible thresholds, the Preisach operator iJ^ is then 
constructed. Memory is represented by means of an antimonotone graph in P, 
which separates relays staying in different states. The reader is advised to refer 
to the appendix for more detailed mathematical description of hysteresis. 

It appears that this mathematical approach provides a realistic approxima- 
tion for the kinetics of ATP and adenosine and their dynamical states, in terms 
of their action on presynaptic terminals. The memory, based on the delayed 
relay operator, characterizes faithfully this biological process. 

2. The Preisach Model 

For the formulation of the ATP kinetics at tripartite synapses in the language of 
hysteresis, the relationship between ATP concentrations and activity behavior 
should be estabilishcd. ATP can be released in concentrations that might be 
subthreshold for the activation of P2X receptors (micromolar) yet sufficiently 
high to allow an accumulation of adenosine that will activate AI receptors (tens 
of nanomolar). 

This means, that the ATP kinetics can be decomposed into two phases: 

1. Hysteron — Ra^M of ATP concentrations acting on P2X receptors, 
ai denotes the releasing point of ATP (concentration zero) and hi is rep- 
resenting the ATP concentration activating P2 purinergic receptors. 
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2. Hystcron = Ra2.b2 of accumulated adenosine (hydrolized ATP) that 
acts inhibitor via Al receptor activation. 02 denotes the maximum con- 
centration level of ATP during pre-accumulation time of adenosine. 62 is 
the accumulated adenosine concentration that is sufficient to activate Al 
receptors. Here, the hydrolization of ATP to adenosine is seen as a change 
in the activity state of ATP; 

The parallel c;oniposition of these hystcrons with equal weights fij leads to the 
finite Weiss-Preisach model of ATP hysteresis: 

2 

y{t) = ^iijR^[to,f]o]x{t) 

where x{t) is the concentration value of astrocytic ATP and rjo the initial 
state of the system. This behavior is shown in figure 1. 

Some important remarks should be mentioned at this point that are suggested 

by experimental data: 

• &2 ~ 10^ • bi. This reflects the concentration thresholds for the activation 

of Al (tens of nanomolar) and P2 receptors (micromolar) . 

• Tf)2 —Tf,^ « 200ms. Ti,. denotes the instant correlating to a certain concen- 
tration. The difference is in accordance with the pre-accumulation time 
of adenosine. 

This simple formulation output variable y{t) of the two-phaseal ATP-induced 
activity provide an useful framework for investigating the purinergic role in 
neuron-glia interactions at tripartite synapses. It also implies that y{t) as a 
hysteresis process, is governed by rate independent memories of ATP concen- 
trations at different instants. 

3. Discussion 

Although this concept is a rough simplification of the ATP-related processes and 
does not provide new insights in the physiology of the tripartite synapses at this 
level, it supplies a powerful tool for modeling of synaptic behavior. It should 
be considered as a pre-step of the modeling the neurochemical interaction in 
tripartite synapses, roughly spoken as a physical characterization that provides a 
useful mathematical relation for further dynamical modeling tasks. The Memory 
of ATP concentrations determines the kinetics of purinergic activity at tripartite 
synapses. The knowledge about this biological memory could be of use for 
physiological investigations. By enlarging the number of hysterons e.g. within 
the knowledge of the release procedure of astrocytic ATP, in terms of vesicular 
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or non- vesicular release, the models refinement and accuracy increases that can 
lead to specific applications of its results other than its usage as a completion 
tool for other models. 

Acknowledgment. . The financial support by the Interdisciplinary Center for Sci- 
entific Computing and the International GraduiertenkoUeg 710 (DFG) and the 
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Appendix. The aim of this appendix, which is merely based on the work of 
Krasnosel'skii and Pokrovskii ([12 ), and Visintin ([19 ), is to introduce the 
general idea of the Preisach model of hysteresis to unexperienced readers. 

The hysteresis nonlinearities are defined as deterministic, rate independent 
operators. The standard approach to hysteresis, is the introduction of ele- 
mentary hysteresis nonlinearities, called hysterons, and treatment of complex 
hysteresis processes as superposition (parallel connection) of these hysterons. 
Among hysterons, the most important are probably the nonideal relays, or ther- 
mostat nonlinearities. 

The nonideal relay (with the threshold values a < 6) is the simplest hysteretic 
transducer. This transducer is denoted by Ra,b- Its output y{t) can take one 
of the two values or 1: at any moment the relay is either 'switched off' or 
'switched on'. The output variable y{t) 

y{t) = Ra.b[h,'no]x{t) , i > io, 

depends on the input variable x{t) e C°([io,T]) and on the initial state 770 € 
{0, 1}. The values of the output function y{t) at an instant t are defined by the 
following explicit formula: 

• ryo, if a < x{t) < bV t e [Lq, t]; 

• 1, if 3 <i e [to,t]: x{ti) > b, x{t) > a V r e [ti,t]; 

• 0, if 3 <i e [to,t]: x{ti) < a, x{t) <by t e [ti,t]; 

Consider a family R of relays i?*" = Ra^,b^ with threshold values aw, b^ (a 
bundle of relays), with w e f2. The index set fl may be finite or infinite. Suppose 
that the set f2 is endowed with a Borel measure /i. Further suppose that the 
functions aw,bw are measurable with respect to fi. We call any measurable 
function ri{w) : il — > {0, 1} the initial state of the bundle R. 

For any initial state rjo^w) and any continuous input x(t), t >to, we define 
the Preisach hysteresis operator as the integral operator: 

y{t) - H^[to,mMt) = [ R'"[to,r]oH]x{t)dfi , t > to . 
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Figure 1: The hysteresis diagramm of the ATP concentration and purinergic receptor activa- 
tion. Two hysterons are parallel connected: First hysteron R^j^^bi is the delay relay operator 
of the accumulated to adenosine hydrolized ATP concentration; The second hysteron Ra2,b2' 
is the delay relay operator of the P2 activating ATP concentrations. The extremal values are 
denoting the concentration threshold of the release point of ATP and the maximum concen- 
tration of tripartite synaptic ATP. 
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